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Abstract. A model with a condensed iron surface and partially ionized hydrogen-thin atmo- 
sphere allows us to fit simultaneously the observed general spectral shape and the broad ab- 
sorption feature (observed at 0.3 keV) in different spin phases of the isolated neutron star RBS 
1223. We constrain some physical properties of the X-ray emitting areas, i.e. the temperatures 
(T poie i ~ 105eV,T poie2 ~ 99 eV), magnetic field strengths (B poi(! i « B po;e 2 ~ 8.6 x 10 13 G) at 
the poles, and their distribution parameters (al ~ 0.61, a2 ~ 0.29, indicating an absence of 
strong toroidal magnetic field component). In addition, we are able to place some constraints 
on the geometry of the emerging X-ray emission and the gravitational redshift (z ~ 0.16^o 3 oi) 
of the isolated neutron star RBS 1223. 
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1. Introduction 

Observations and modeling of thermal emission from isolated neutron stars can provide 
not only information on the physical properties such as the magnetic field, temperature, 
and chemical composition of the regions where this radiation is produced, but also we 
may infer on the properties of matter at higher densities deeper inside the star. 

In particular, the study of thermal emission from isolated neutron stars may al- 
low one to infer the surface temperature and total flux measured by a distant ob- 
server and to estimate the real parameters. With the known distance and the red- 
shifted radius of the neutron star the actual radius and mass of a neutron star are: 
R = R°°[l-2GM/Rc?} 1 / 2 , M= ^§ [l - (-J^) 2 (see, e.g. Zavlin (2009)) 

The detection and identification of any absorption/emission feature in the spectrum or 
performing rotational phase-resolved spectroscopy of isolated neutron stars will allow us 
to determine gravitational redshift and directly estimate the mass-to-radius ratio, M/R. 
Together they yield a unique solution for M and R. These spectral features may allow 
to measure the neutron star magnetic field and provide an important input for modeling 
of magnetized atmospheres. 

RBS 1223 shows the highest pulsed fraction (13-42%, depending on energy band, see 
Fig. [1]) and strongest broad absorption feature (Schwope et al. 2007) of all isolated 
neutron stars. 

High quality rotation phase resolved spectroscopy is needed in order to fit with neutron 
star magnetized atmosphere models and to constrain the gravitational redshift of the 
neutron star (Suleimanov et al. 2010; Hambaryan et al. 2011). 
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Figure 1. XMM-Newton EPIC pn co-added phase-averaged X-ray spectra (left panel) including 
primary and secondary peaks, first and second minima, and phase-folded light curves (right 
panel) in different energy bands of RBS f 223 combined from f 2 pointed observations. 



2. Data analysis and results 

Using the data collected with XMM-Newton EPIC pn from the 12 publicly (similar in- 
strumental setup, i.e. Full Frame, Thinl Filter) available observations, in total presenting 
about 175 ks of effective exposure time, we extracted spin-phase resolved spectra with 
high S /N ratio (Fig. [IJ and fitted simultaneously with highly magnetized isolated neu- 
tron star surface/atmosphere models (Suleimanov et al. 2010). These models arc based on 
various local models and compute rotational phase dependent integral emergent spectra 
of isolated neutron star, using analytical approximations. 

The basic model includes temperature/magnetic field distributions over isolated neu- 
tron star surface^, viewing geometry and gravitational redshift. Three local radiating 
surface models are also considered, namely, a naked condensed iron surface (van Adels- 
berg & Lai 2006) and partially ionized hydrogen model atmospheres, semi-infinite or 
finite atop of the iron condensed surface. 

The observed phase resolved spectra (i.e. energy-spin phase image, Fig. [2]) of the iso- 
lated neutron star RBS 1223 are satisfactorily fitted (verified also via MCMC, Fig. [3]) 
with two slightly different physical and geometrical characteristics of emitting areas, a 
model parameterized with a Gaussian absorption line superimposed on a blackbody spec- 
trum and by the model of a condensed iron surface, with partially ionized, optically thin 
hydrogen atmosphere above it, including vacuum polarization effects, as orthogonal ro- 
tator. Note, the latter one is more physically motivated. We have additionally performed 
Markov Chain Monte Carlo (MCMC) fitting as implemented in XSPEC. 

The fit also suggests the absence of a strong toroidal magnetic field component. More- 
over, the determined mass-radius ratio, (M/Ms„„)/(i?/km) = 0.087 ± 0.004, suggests a 
very stiff equation of state of RBS 1223. 

Acknowledgement 

f T 4 = Tpi,2 cos a I sin a e + T tm, B = B p i, 2 \J cos 2 9 + a ia sin 2 9, where the parameters a l:2 
axe approximately equal to the squared ratio of the magnetic field strength at the equator to the 
field strength at the pole, ai,2 ~ (B C q/ -B P i,2) 2 • Using these parameters we can describe various 
temperature distributions, from strongly peaked (a ^ 1) to the classical dipolar (a = 1/4) and 
homogeneous (a = 0) ones. 
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Figure 2. Energy- Phase image 
of RBS1223 combined from dif- 
ferent observations. Rotational 
phase-folded light curve in the 
broad energy, 0.2-2.0 keV, band 
(left panel) and the phase aver- 
aged spectrum (bottom panel) 
are shown. 




Figure 3. Probability density distribution of gravitational redshift by Markov Chain Monte 
Carlo (MCMC) fitting with the model of strongly magnetized neutron star condensed Iron 
surface and partially ionized Hydrogen thin atmosphere atop it (left panel). Dashed vertical 
lines indicate the highest probabil ity interval (68%). Mass-radius relations (right panel) for 



several EOS (Haensel et al. (2007) thin solid curves), and a strange star (thin dash-dotted). 



Thick dashed: curve of constant 7?°° = R/y/l - 2GM/Rc 2 = 17 km ( |Trumper (20 05) 
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